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Abstract. The influence of antiferromagnetic order on the mixed state of a superconductor may result in 
creation of spin-flop domains along vortices. This may happen when an external magnetic field is strong 
enough to flip over magnetic moments in the vortex core from their ground state configuration. The 
formation of domain structure causes modification of the surface energy barrier, and creation of the new 
state in which magnetic flux density is independent of the applied field. The modified surface energy barrier 
has been calculated for parameters of the antiferromagnetic superconductor DyMogSs. The prediction of 
two-step flux penetration process has been verified by precise magnetization measurements performed on 
the single crystal of DyMogSs at milikelvin temperatures. A characteristic plateau on the virgin curve 
B(Ho) has been found and attributed to the modified surface energy barrier. The end of the plateau 
determines the critical field, which we call the second critical field for flux penetration. 

PACS. 74.60.-w Type-II superconductivity - 74. 25. Ha Magnetic properties 



Introduction 

The discoveries of ternary Rare Earth (RE) Chevrel Phases 
REMogSs and RERI14B4 compounds with regular dis- 
tribution of localized magnetic moments of RE atoms have 
proved conclusively the coexistence of various types of 
magnetism with superconductivity. Intensive experimen- 
tal and theoretical works have shown that 4f electrons 
of RE atoms responsible for magnetism and 4d electrons 
of molybdenum chalcogenide or rhodium boride clusters 
responsible for superconductivity are spatially separated 
and therefore their interaction is weak. It seems that in 
many of these systems superconductivity coexists rather 
easily with antiferromagnetic order, where usually the Neel 
temperature Tjv is lower than the critical temperature for 
superconductivity T c . On the other hand, ferromagnetism 
and superconductivity cannot coexist in bulk samples with 
realistic parameters. Quite often the ferromagnetic order 
is transformed into a spiral or domain-like structure, de- 
pending on the type and strength of magnetic anisotropy 
in the system ||,||. For almost two decades the problem 
of the interaction between magnetism and superconduc- 
tivity has been overshadowed by high temperature super- 
conductivity found in copper oxides. However, the recent 
discovery of the presence of magnetic order in Ru-based 
superconductors has triggered a new series of ex- 

periments and inspired a return to the so-called coexis- 
tence phenomenon Q. Most recently, the interplay be- 
tween magnetism and superconductivity was studied in 
d-electron UGe2 J§| and ZrZn 2 |J, where itinerant fer- 
romagnetism may coexist with superconductivity, and in 



heavy fermion UPd2Al3 [[To), where magnetic excitons are 
present in superconducting phase. 

Among classical magnetic superconductors, the Chevrel 
phases have been studied most intensively. These com- 
pounds are mainly polycrystalline materials. However, some 
specific features can be measured only on single crystals. 
One such effect is a two-step flux penetration process, pre- 
dicted in Ref. ( [ pT|p^ ] ) and later observed in an antiferro- 
magnetic superconductor (bet) ErRli4B4 [[l3| and recently 
in DyMogSs fTi)] . This anomaly was explained as a result 
of the magnetic transition taking place in the vortex core. 
This transition seems to create a new type of vortices with 
the unique magnetic structure as shown in FigQ In the 
present paper the two-step flux penetration process is cal- 
culated and measured for a single crystal of DyMogSs- 
The DyMoeSs compound with T c = 1.6 K exhibits transi- 
tion from the paramagnetic to the antiferromagnetic state 
at Tn = 0.4 K. Its crystal structure can be described as 
interconnected Mo 6 Ss units and Dy ions. One such unit 
is a slightly deformed cube where S atoms sit at the cor- 
ners and Mo atoms are situated at the cube-faces. The 
MogSs units are arranged in a simple rhombohedral lat- 
tice and Dy ions are located in the center of the unit cell. 
The magnetic moments of Dy ions form a simple structure 
consisting of (100) planes with moments of 8.7 hb alter- 
nately parallel and antiparallel to the [111] rhombohedral 
axis. 

Neutron experiments performed on DyMogSs in an ap- 
plied magnetic field at T = 0.2 K have revealed in the 
intensity spectrum a number of peaks characteristic for 
ferromagnetic order jl5). These peaks begin to develop 
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Fig. 1. The magnetic structure of the isolated vortex and the 
distribution of the magnetic induction around the vortex core 
in the spin-flop SF and antiferromagnetic AF phases jnj]. 



at Hq — 200 Oe, much below the superconducting upper 
critical field H C 2- Thus, in DyMogSs a kind of ferromag- 
netic order coexists with superconductivity in the same 
manner as antifcrromagnetism. For a field applied paral- 
lel to the [111] direction (magnetic easy-axis direction), 
the ferromagnetic order is a spin- flop type ]16[ . This fea- 
ture is easy to understand. Consider the well known phase 
diagram of a two-sublattice antiferromagnet. An infinites- 
imal magnetic field applied perpendicular to the easy axis 
makes the ground antiferromagnetic configuration unsta- 
ble against the phase transformation to the canted phase. 
On the contrary, if the magnetic field is applied parallel 
to the easy axis the antiferromagnetic (AF) phase is sta- 
ble up to the thermodynamic critical field Ht- When the 
field is further increased, a spin-flop (SF) phase develops 
in the system. Let us assume that in the antiferromagnetic 
superconductor the lower critical field fulfills the relation 
H c i < \Ht and that the external field Hq is applied par- 
allel to the easy axis. When H cl < Hq < \H T the vortices 
appear entirely in the AF phase. When Hq is increased be- 
yond \Ht the phase transition to the SF phase originates 
in the core, because near H r j t he field intensity in the core 
is approximately twice H c \ |17[ . The spatial distribution of 
the field across the vortex is a function decreasing from the 
center as seen in Fig.|l|. Thus, the magnetic field intensity 
outside the core is less then Ht and, therefore, the rest of 
the vortex remains in the AF phase. The radius of a spin- 
flop domain grows as the field is increased. The formation 
of domains inside the vortices should be accompanied by 
the modification of the surface energy barrier |12|] . This 
process leads to a state of the superconductor in which 
flux entrance is temporarily prohibited - flux density is 
independent of the applied field. In order to kill this state 
the external field should be increased above certain second 
critical field for flux penetration. Then, the vortices pene- 
trating the sample will have the spin flop domains created 
along the cores. The goal of our work is to compare the 
model calculations based on the method of images |Q8| 



with the experimental results obtained for the DyMo@Sg 
single crystal. 



Theoretical considerations 

In order to describe thermodynamic behavior of DyMogSs 
for constant temperature and constant external magnetic 
field we introduce the following free energy functional 



F= I dV{fs + f M + ^-(B-4nM) 2 }. 

J 87T 



(1) 



Here B is the magnetic induction, M = Mi + M 2 the 
magnetization of two sublattices antiferromagnet and H = 
B — 47rM the intensity of thermodynamic magnetic field. 
The energy density of the superconducting subsystem fs 
is expressed in a standard way 



2m 



cn 



am 



(2) 



The quantity e,m,c denote charge and mass of the elec- 
tron and light velocity, respectively. According to experi- 
ments the antiferromagnetic order is practically unaffected 
by the appearance of superconductivity, then it is rea- 
sonable to neglect the effect of superconductivity on the 
exchange interaction in F. This means that both order 
parameters <P and M are coupled via the vector potential 
A: 



VxA = B = H + 4ttM 
j s = -£-V x H, 



(3) 



where j s denotes a superconducting current. The antifer- 
romagnetic energy density, which mimics the experimental 
results in DyMogSg, is given by the following expression 



2 



f M = JM 1 -M 2 +K^(M?) 

i 



2 



i=l j=x,y,z 



(4) 

Here J is the exchange constant between two antiferro- 
magnetic sublattices, K denotes single ion anisotropy con- 
stant, and y^pyf is the magnetic stiffness length. Since in 
the following we analyze the phenomena with character- 
istic length-scales much larger then the interatomic Dy- 
Dy distance it is justified to omit the gradient term in 
/m- The components of the total magnetization vector 
M = Mi + M 2 , | Mi | = M (i = 1,2) have the follow- 
ing form in both sublattices: Mix = Mo sin ^i, Mi y = 
0, M iz = M cos$i, where 9i (canted spin angle) is the 
angle between the magnetization in the sublattice and 
the external magnetic field directed along z-axis. The AF 
(#i = 0, 8 2 = tt) and SF phases (9i = -9 2 = 8) are in 
thermodynamic equilibrium in an applied field equal to 
the thermodynamic critical field. 



H T = M [K(J - K)] 



1/2 



(5) 
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The canted spin angle of the SF phase is then expressed 
as 

cost* = — — . (6) 
tlx 

The equilibrium conditions of the whole system can be ob- 
tained via minimization the Gibbs free energy functional 

G = F — f(B ■ H W with respect to <f , A and M. 

4?r J. 

Performing this task in London approximation one can 
obtain 

B + A 2 V x V x (B — 4ttM) = 0, (7) 

where A is the London penetration depth. The appropri- 
ate equations describing spatial distribution of M should 
accompany Eq.(Q). To make the problem simpler we sup- 
pose that the magnetization is constant across the SF do- 
main [|l2],[n| ■ In this way the distribution of the magneti- 
zation around a single vortex is the following 



a 



1o 
2ttA 2 



+ tMtMt)P 



(12) 



Finally, the minimum of the free energy of the vortex per 
unit length 



£"1 



A 2 

8tt 



f dl {[b S F - 4?rM] x V x b S F} 
<b dl {b A F x V x b A F} , 

J <Jo 



with respect to rg determines: 



IMI 



M if r < ro 
if r > r ' 



(8) 



ro _ 

A J ~ ttX 2 B t ' 



where ro is the radius of the spin-flop domain. With the 
help of Eq. (Q) one can write Eq. (0) for a single vortex as 
follows 

F = 8^ / {^ SF ~ 4?rM ) 2 + a2 I v x ( b SF - 4ttM)] 2 | dV SF 



8tt 
1 

8^ 



b| F + A 2 (V x b AF r 



dV AF . 



(9) 



Here b A p and bgp denote magnetic induction in AF and 
SF phases of a single vortex, respectively. The integrals are 
performed over the volume of each phase with the exclu- 
sion of the volume of the vortex core. Equation (Q), for a 
single vortex, can be solved in the cylindrical coordinates 
in terms of the modified Bessel functions K$ and Iq : 

b SF = dK (j) + C 2 I Q (£) , for £ < r < r 
b AF = C 3 K (£) , forr> r , (10) 



(£ denotes the coherence length) with the following bound- 
ary conditions: 



6 S f (y ) =H T + 4ttM = B T 

&AF (y) = 



(13) 



(14) 



The line integrals in Eq. ( pL3| ) are performed over the cross- 
sections perpendicular to the axis of an appropriate cylin- 
drical element of the surface of the vortex, <j\ denotes the 
surface of the vortex core, a 2 the surface of the SF domain. 

In order to study the conditions under which magnetic 
flux density in the sample becomes unstable in the applied 
magnetic field one must take into account the surface en- 
ergy barrier preventing vortices from entering or exiting 
the sample. The presence of a surface of the supercon- 
ductor leads to the distortion of the field and current of 
any vortex located within a distance of the order of pene- 
tration depth from the surface. The requirement that the 
currents cannot flow across the surface of the supercon- 
ductor leads to the introduction of an image vortex, at 
x = —xl, with vorticity opposite to the real one. Both 
vortices interact as real ones except that the interaction 
is attractive. 

We consider semi-infinite specimen in the half space 
x > 0, the vortex and the external magnetic field run- 
ning parallel to the surface. In the low flux density regime 
£ 2 < (fo/B < A 2 , Clem |Q has shown the existence of a 
vortex-free region of the width x v f near the surface of the 
sample and constant vortex density region for x > x v f. 
Within vortex-free area one can introduce locally aver- 
aged magnetic field Bm exponentially decreasing from its 
surface value H to its average interior value B, 



(11) 



B 



M 



B cosh 



A 



(15) 



These conditions, together with the flux quantization con- 
dition, are used to calculate the arbitrary constants in 
Eq.©. 



The boundary condition Bm(0) — Hq determines the thick- 
ness of the vortex-free region 



Ct = 

c 2 = 
c 3 = 







B T -h 



(?) 



- alo 



X v f 



A cosh ' ( —p 
B 



(16) 
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Now we can characterize the distribution of the magnetic 
induction around a single vortex in the vortex-free region 



■•(?) 
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(17) 



The Gibbs free energy of the system can now be written 
in the following way 

A 2 f 

G = — <h der { [Bsf - 2H - 4ttM] x V x B S f} 

87r J<T1 

A 2 r 

+ — A der {[Baf - 2H ] x V x B AF } 
A 2 f 

+ —j der {zB M x V x B AF } , (18) 

where 173 is the surface of the specimen. After some trans- 
formations Jl2|,[l8), one can obtain the Gibbs free energy 
per unit length Q in the following form: 



A 2 
2 



D!H - D 2 B M 



X y J X 

A 



(19) 



where 



db 



SF 



d x = - e- 



x - X L 



dx 



dbsi 



ro- 



x - x L 



dx 



db af 



ni- 



di) 



X - X L 



SF 



Do 



dx 



x - x L 



dx 



db S t 



ro- 



x- x L 



dx 



C?&AF 



2r - 



X - XL 



dx 



x' = x L + £ ; x" — x L + r 

Q has its maximum at x = x max somewhere in the vortex- 
free region ro < x max < x v f . We can find x max solving the 
force balance equation. When the external field reaches 



H en2 (B) =5 cosh (^p 



(20) 



where x en is the vortex-free width corresponding to an ex- 
ternal field equal to H en 2 , the energy barrier moves within 



a distance tq of the surface (ro <C x v f). Thus, one can get 



XDi 
2Lh 



db 



AF 



dx L 



= B sinh 



'•() 



(21) 



The left hand side of the above equation gives H en 2(0). 
This field may be thought as the second critical field for 
flux penetration calculated in the single vortex approxi- 
mation fill. 



2H en2 {0) = 



H 7 



ipo , / it\ 2 Bt 



■In 



(22) 



it\ 2 Bt \ <po 
Taking into account that ro <ti x en we finally obtain 



H en2 (B) = JB 



(23) 



Let us make a short summary of the calculations. When 
the SF domain develops, the screening current must redis- 
tribute its flow around the vortex in order to fulfill the sin- 
gle flux quantum requirement. This one can easy deduce 
from Eqs (10)-(p~2|). The redistribution of the screening 
current causes the change in the surface energy barrier 
preventing vortices from entering into the sample. This 
is expressed in Eq.(19). Consequently, the averaged flux 
density in the sample B = nip® remains constant when 
the external field is increased. The vortices start to pene- 
trate into the sample again when the second critical field 
for flux penetration, calculated in Eq.pSj), is reached. 



Experimental details 

The single crystals of DyMogSs were grown using the 
slow cooling of a melted charge closed in hermetically 
sealed molybdenum ampoules. Details of the crystal grow 
procedure are described elsewhere |l^,^0|. The crystals 
were pure, homogeneous and large enough to be used for 
studying some subtle effects accompanying the magneti- 
zation process at milikelvin temperatures. Chemical com- 
position and crystal uniformity were examined using a Hi- 
tachi Scanning Electron Microscope equipped with an en- 
ergy dispersive x-ray analyzer. Single-crystal x-ray diffrac- 
tion measurements were performed at room temperature 
on a Simens SMART CCD diffractometer. The electron 
probe microanalysis of the regular-shaped crystals showed 
a composition corresponding to the DyMoeSs formula. 
The cell parameters in the rhombohcdral lattice were a^t = 
6.452 • 1CP 8 cm and an — 89.50°, and were equivalent for 
all crystals analyzed. The single crystal selected for our 
experiment had dimensions 0.2 x 0.2 x 0.2 mm 3 and a 
mass ~ 0.05 mg. 

Magnetization was measured with the SHE 330X se- 
ries SQUID system with SQUID sensor installed in the 
vacuum chamber of the 3 He - 4 He dilution refrigerator. 
The sensor was thermally anchored to the liquid He bath 
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(4.2 K) and shielded with a Nb tube. Two counter-wound 
pickup coils were connected to the input coil of the SQUID 
sensor. The SQUID pickup coils were placed in the cen- 
ter of a 10 cm long superconducting solenoid generating 
a magnetic field up to 1.5 kOe. Both the coils and the 
solenoid were fixed to the mixing chamber of the dilu- 
tion refrigerator. Details of the experimental setup are de- 
scribed elsewhere [14]. The perfect shielding (iwAi = Hq) 
of the Meissner state was used to calibrate the SQUID 
system. The crystal was oriented with the magnetic easy 
axis (the [111] crystallographic triple axis) parallel to the 
external magnetic field. For this orientation, the demag- 
netizing factor was assumed to be k = 1/3. 



Comparison with theory 

In Fig.^J, the magnetization M. measured as a function 
of temperature is presented for several applied magnetic 
fields oriented parallel to the easy axis of the single crys- 
tal. At higher temperatures, the transition to the super- 



•*—> 
'c 

is 
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Fig. 2. Magnetization versus temperature at several applied 
fields for DyMogSs single crystal for the field direction oriented 
parallel to the magnetic easy axis. 



conducting state is observed at T c as the smooth decrease 
of M (e.g., T c = 1.62 K for H = 20 Oe). This critical 
temperature is clearly field dependent as expected for a 
superconductor. At low fields, M. riches negative values 
close to T c . At higher fields, this is not possible because of 
the induced strong paramagnetic moment of the Dy ions. 
At low temperatures, the abrupt change of M. is observed 
at Tjv = 0.4 K, reflecting the transition to the AF state. 
In that state, the internal field is reduced and M. can 
now become negative even for higher fields. At Tjv and for 
Hq < 200 Oe, the change of M. between the paramagnetic 
(PM) and AF states increases significantly with increasing 



field, as expected. However, for H > 200 Oe, the single 
crystal is in the SF phase Jl3,|l6| and the observed change 
of Ai, caused by the transition to the ordered state, now 
decreases with increasing field. 



200 




300 



H (Oe) 





Fig. 3. Magnetic induction for DyMogSs single crystal in the 
virgin state measured as a function of an applied field for three 
temperatures below Tjv = 0.4 K. The field direction is oriented 
parallel to the magnetic easy axis of the crystal. Each B(Hq) 
curve exhibits characteristic plateau indicating that a number 
of vortices is kept constant when the external field is increased. 
The results are not corrected for demagnetizing effects. The 
corrected values used for calculations are given in Table ^. The 
solid lines are guide to the eye. 



The low- field parts of the B(Hq) virgin curves are pre- 
sented in Fig.|| to show the details of flux penetration. 
These curves have been obtained by the simple trans- 
formation of the M(Hq) results reported in our previ- 
ous work 1 14 . The observed penetration is typical above 
Tjv = 0.4 K and proceeds as an unusual two-stage process 
at lower temperatures where AF order coexists with su- 
perconductivity. At low fields the sample is in the Meissner 
state. When the field increases above H c i, the sample is 
penetrated by the flux. Then, at higher fields, the penetra- 
tion process stops unexpectedly and B = B p i is constant 
in the sample while the external field is further increased. 
This new perfect shielding appears at H = H p i. The 
penetration process starts again when the field reaches 
H = H en 2- This value we call the second critical field for 
flux penetration. 

An interesting effect of the temperature dependence of 
H c i below is seen in Fig. 3 and Table 1. The very small 
decrease in temperature results in the significant increase 
of H c \, whereas i? C 2 measured by us in 0.1 K, 0.12 K, and 
0.14 K remains close to 900 Oe. The similar behavior has 
been observed in GdMogSes, its main feature is a sharp 
dip on the H c i(T) around TV and a plateau on H C 2(T) 
for T < TV f21|] . This behavior agrees well with our 
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theoretical model. The large increase of H c \ for tempera- 
tures decreasing below TV indicates that our sample is in 
the antiferromagnetic collinear phase. In this phase, when 
temperature decreases from Tjv, the pair-breaking effects 
due to molecular field and antiferromagnetic fluctuations 
weaken and H c i(T) rapidly tends to match its pattern in 
the paramagnetic state. On the contrary, the nearly con- 
stant H C 2 below Tjv |^2|,|2^] indicates that the sample is in 
the antiferromagnetic canted phase and pair-breaking ef- 
fects due to the on-field component of the molecular field 
are present even at the lowest temperatures. In our model 
we have assumed that vortices enter the sample in the 
collinear antiferromagnetic state, and when the field is in- 
creased the canted phase appears first inside the core of 
the vortex. Thus, above H en 2, if the external field is fur- 
ther increased the volume of the canted phase enlarges. 
This makes possible to transform the whole sample to the 
canted phase well below H c2 . 

In order to compare our theoretical model with the 
experimental results we have estimated several quantities. 
The most important is the Ginzburg-Landau (GL) param- 
eter k. This parameter has been calculated in p4| to be 
equal to 2.6 for T — 0.10 K. We have taken advantage of 
the strange behavior of H c i(T) and H C 2{T) to calculate 
k for T = 0.14 K and T = 0.12 K. The constant value of 
H c2 = 900 Oe in the range of 0.10 K < T < 0.14 K pre- 
dicts that coherence length does not change in this interval 
of temperatures where an abrupt increase of H c i suggests 
that the penetration depth drastically decreases. This ob- 
servation has been used to write the following equation 



H cl (T ) = (k{T)V 1h«(Tq) 
H cl (T) U(W ln«(T)' 



(24) 



where T = 0.10 K , k(T ) = 2.6 and 0.12 K < T < 
0.14 K. The above equation has been solved numerically 
and the results for k are given in the Table 



Table 1. Experimental values of Fig.^| corrected for demagne- 
tizing effects according to the formulae: H — Ho + 4nkM, and 
B — H — 47r(l — k )M p4| ], where M (absolute value) is taken 
from Fig.4 of Ref.O. 



To find the thermodynamic critical field Ht and then 
to calculate H en 2(B) we have used the following argumen- 
tation. At low fields, in the vicinity of the lower critical 
field H c i, the intensity of the field in the vortex core is 
2iJ c i Jl7j . When the external field is increased the field 
intensity in the vortex core increases because of the su- 
perposition of the fields of the surrounding vortices. The 
field intensity in the core must reach Ht in order to orig- 
inate the transition to the SF phase. Thus, taking into 



account only the nearest z neighbors we can write 



Ht = 2H c1+Z ^-kJ\ 
2irX z \ X 



(25) 



where d denotes intervortex spacing, and d/X corresponds 
to the value B p i for which the penetration process unex- 
pectedly stops. The relations Ba = 2tp /d 2 V3 (for tri- 
angular lattice of vortices), tp = 2ttH C 2^ 2 , where H C 2 = 
900 Oe [[l4| have been used to obtain d/X. Then, this 
value has been inserted into Eq.(^5|) to obtain Ht- The 
saturation magnetization of Dy ions, 8iyMq = 3780 G has 
been calculated taking into account the volume of the el- 
ementary cell of DyMo 6 S 8 , V = 268 • 10~ 24 cm 3 . The 
anisotropy coefficient K = 0.44 has been determined for 
each magnetization curve by finding the best fit of the the- 
oretical with the experimental magnetization curves p5[ . 
Next, the magnetization in the SF-phase domain has been 
calculated with the help of Eq. (||) 



M = 2M o cos0 



2KM 2 
H T 



(26) 



Eq.(p6|) gives M corresponding to the field H p \ for which 
the penetration stops. Finally, inserting all the above cal- 
culated values into Eqs.( [22|]23| ) we have obtained H en2 (B). 
The results are summarized in the Table @. 



Table 2. Summary of the calculated quantities. 



T [K] 


d/X 


Ht [Oe] 


4tvM [G] 


B T [G] 


H m2 (B) [Oe] 


0.14 


1.7 


250 


1000 


1250 


215 


0.12 


2.8 


325 


775 


1100 


240 


0.10 


4.9 


360 


700 


1060 


265 



T [K] 




Hd [Oe] 


H pl [Oe] 


B pl [G] 


H en2 (B) [Oe] 


0.14 


4.3 


100 


170 


135 


250 


0.12 


3.1 


150 


185 


105 


270 


0.10 


2.6 


180 


200 


80 


280 



Conclusion 

We have demonstrated that the antiferromagnetic super- 
conductor DyMogSs shows interesting behavior in the mag- 
netic field applied below T^. The sample in the virgin 
state magnetizes initially like ordinary type II supercon- 
ductor. When the applied magnetic field reaches the crit- 
ical field for flux penetration the sample transforms from 
the Meissner to the mixed state. Then, magnetization pro- 
ceeds in an unusual way. As the field is further increased, 
a new shielding state appears but, in the contrary to the 
Meissner state, with a constant flux density inside the 
sample. Characteristic plateau, observed for the magne- 
tization curves, proofs that magnetic flux density inside 
DyMo 6 Ss is unaffected by the increased external field. 
When the field reaches certain value, we call it the second 
critical field for flux penetration, the flux starts to enter 
the sample again. This phenomenon we name two-step 
flux penetration. We have argued that in this new state 
vortices transform to the shape shown in Fig.^, where a 
domain of the spin-flop phase is created. The expected 
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metamorphosis of the vortices leads to a spatial redistri- 
bution of the shielding supercurrents, flowing around the 
core, in order to keep constant the flux carried by each 
vortex. Consequently, a new energy barrier is formed near 
the surface preventing vortices from entering the sample. 
Thus, the number of vortices inside the superconductor is 
kept constant. To overcome the new energy barrier by the 
vortices with magnetic structure the external field must be 
increased beyond H en2 , the second critical field for flux 
penetration. The formula for this field has been derived 
using the image method. The values of H en2 calculated 
for three temperatures below Tm agree very well with the 
experimental results. 
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